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The nature of the glass transition, whereby a variety of liquids transform into 
amorphous solids at low temperatures, is a subject of intense research despite 
decades of investigation. Explaining the enormous increase in viscosity and re- 
laxation times of a liquid upon supercooling is essential for an understanding 
of the glass transition. Although the notion of a growing length scale of "coop- 
eratively rearranging regions" (CRR), embodied in the celebrated Adam-Gibbs 
relation, has been invoked for a long time to explain dynamical slow down, the 
role of length scales relevant to glassy dynamics is not well established. Recent 
studies of spatial heterogeneity in dynamics provide fresh impetus in this direc- 
tion. A powerful approach to extract length scales in critical phenomena is finite 
size scaling, wherein a system is studied for a range of sizes traversing the length 
scales of interest. We perform finite size scaling for the first time for a realistic 
glass former, using computer simulations, to evaluate a length scale associated 
with dynamical heterogeneity which grows as temperature decreases. However, 
relaxation times which also grow with decreasing temperature, do not show the 
same kind of scaling behavior with system size as the dynamical heterogeneity, 
indicating that relaxation times are not solely determined by the length scale 
of dynamical heterogeneity. We show that relaxation times are instead deter- 
mined, for all studied system sizes and temperatures, by configurational entropy, 
in accordance with the Adam-Gibbs relation, but in disagreement with theoret- 
ical expectations based on spin-glass models that configurational entropy is not 
relevant at temperatures substantially above the critical temperature of mode 
coupling theory. The temperature dependence of the heterogeneity lengthscale 
shows significant deviations from theoretical expectations, and the length scale 
one may extract from the system size dependence of the configurational en- 
tropy has much weaker temperature dependence compared to the heterogeneity 
length scale at all studied temperatures. Our results provide new insights into 
the dynamics of glass-forming liquids and pose serious challenges to existing 
theoretical descriptions. 

Most approaches to understanding the glass transition and slow dynamics in glass form- 

ers 0, y, y, y, y, y, y, y, 3 are based on the intuitive picture that the movement of their 

constituent particles (atoms, molecules, polymers) requires progressively more cooperative 
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rearrangement of groups of particles as temperature decreases (or density increases). Struc- 
tural relaxation becomes slow because the concerted motion of many particles is infrequent. 
Intuitively, the size of such "cooperatively rearranging regions" (CRR) is expected to in- 
crease with decreasing temperature. Thus, the above picture naturally involves the notion 
of a growing length scale, albeit implicitly in most descriptions. The notion of such a length 
scale, related to the configurational entropy S c (see methods), forms the basis of rationaliz- 
ing IlL the celebrated Adam-Gibbs (AG) relation |l| between the relaxation time and 

s c . 

More recently, a number of theoretical approaches have explored the relevance of a grow- 
ing length scale to dynamical slow down . A specific motivation for some of these 



approaches arises from the study of hetero gene ous dynamics in glass formers nj, 11, 12, 13]. 
In particular, computer simulation studies [ill. Il2l . Il3| ] focused attention on spatially corre- 
lated groups of particles which exhibit enhanced mobility, and whose spatial extent grows 
upon decreasing temperature. The spatial correlations of local relaxation permits identi- 
fication of a dynamical (time dependent) length scale, £, through analysis of a four point 
correlation function first introduced by Dasgupta et al 14J (see methods), and the associated 



dynamical susceptibility Xa 



16] . These quantities have been studied recently via inhomo- 



data 
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geneous mode cou plin g theory (IMCT) |5|] and estimated from simulation and experimental 
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The method of finite size scaling, used extensively in numerical studies of critical phe- 
nomena [22 j . is uniquely suited for investigations of the presence of a dominant length scale. 
This method involves a study of the dependence of the properties of a finite system on its 
size. We study a binary mixture of particles interacting via the Lennard- Jones potential 23] ]. 
originally proposed as a model for iVigo-f^cb an d widely studied as a model glass former. We 
perform constant temperature molecular dynamics simulations at a constant volume (see 
methods and 2jj] for details), for seven temperatures, and up to a dozen different system 
sizes for each temperature. For each case, we calculate the dynamic susceptibility Xi{t) as 
the second moment of the distribution of a correlation function Q(t) which measures the 
overlap of the configuration of particles at a given time with the configuration after a time t 
(see methods). From previous work, it is now well-established that Xa{^) has non-monotonic 
time dependence, and peaks at a time T4 that is proportional to the structural relaxation 
time r. Such behavior is shown in the inset of Figure 1(a). In Figure 1(a), we show the peak 



3 



values x\ = X4( T i) vs system size (number of particles) A for a range of temperatures. At 
each temperature, x\ is an increasing function of A, saturating at large A. The saturation 
occurs at a larger value of A at lower temperatures. This is precisely the finite-size scaling 
behaviour expected for a quantity whose growth with decreasing temperature is governed 
by a dominant correlation length that increases with decreasing temperature. 

We have estimated the correlation length £ from finite size scaling of Xa(T,N), which 
also involves estimating the value of x\ as N — > 00. As the latter estimation is a potential 
source of error in estimating £, we employ the Binder cumulant of the distribution of Q(t$) 
to estimate £. The Binder cumulant 25J, defined (see methods) in terms of the fourth and 
second moments of the distribution, vanishes for a Gaussian distribution, whereas it acquires 
negative values which are larger for distributions that deviate from a Gaussian more strongly. 
The Binder cumulant has been used extensively in finite size scaling analysis in the context 
of critical phenomena, owing to its very useful property that in systems with a dominant 
correlation length £, it is a scaling function only of L/£ (or equivalently, of A/£ 3 ), where 
L is the linear dimension of the system. The distributions themselves are shown in the 
inset of Figure 1(b), for two different system sizes for temperature T = 0.47. We see that 
the distribution is unimodal for the large system size of A = 1600 whereas it is strongly 
bimodal for the small system size of A = 150. The same trend is observed as temperature 
is decreased for a fixed size of the system. The data collapse of the Binder cumulant, 
from which we extract the correlation length £(T), is shown in Figure 1(b). The collapse 
observed is excellent, confirming that the growth of x\ with decreasing T is governed by a 
growing dynamical correlation length. The values of £ obtained from this scaling analysis are 
consistent with less accurate estimates obtained from a similar analysis of the A-dependence 
of x\{T, A), as well as from the wave vector dependence of the four point dynamic structure 
factor Si(q,Ti) (see, e. g. [5j). Estimated x\ as A — > 00 compare very well with the 
q — > limit of S±(q, T4), up to a proportionality constant, implying that concerns about the 
dependence of x\ values on the details of the simulations do not apply in our case [26]. 

The value of £ (up to a temperature-independent multiplicative factor that can not be 
determined from scaling, but has been fixed to match results obtained from S^(q,T^)) in- 
creases from 2.11 to 6.23 as T decreases from T = 0.70 to T = 0.45. We find that both £ 
and the asymptotic, A — >• 00 value of x\ deviate from power law behavior as the critical 
temperature T c of mode coupling theory (T c ~ 0.435 in our system) is approached (consis- 
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tently with previous observations). However, the power-law relationship between x\ an d £ 
predicted in IMCT is approximately satisfied by our data. Since the range of the measured 
values of £ is small, it is not possible to obtain accurate estimates of the exponents of these 
power laws. 

Next we consider the dependence of the relaxation time r on T and N. For each case, we 
calculate the relaxation time from the decay of < Q(t) >. The results for r are displayed in 
Figure 2, which shows that r increases as the temperature decreases, as expected. However, 
the observed increase in r with decreasing N for small values of N at fixed T is not consistent 
with standard dynamical scaling for a system with a dominant correlation length (e.g. near 
a critical point): dynamical finite-size scaling would predict a decrease in r as the linear 
dimension L of the system is decreased below the correlation length £. Similar finite size 
effects on relaxation times have been observed in previous simulations of realistic glass 
formers (e. g. j^]) but have not been analyzed in detail. The observed ^-dependence of 
r is opposite to that found in finite-size scaling studies of mean-field spin-glass models 



2a 

291 ] whose analytically tractable behaviour forms the starting point of most recent theories 
(including IMCT) of the dynamics of glassy liquids 30]. 

The inset of Figure 3 shows the large- iV value of r plotted as a function of (a) the 
correlation length £ on a double-log scale, and (b) jA^ on a semi-log scale. The power- 
law relation between these two quantities predicted in IMCT j^J is found above T = 0.5; 
deviation from a power law is found at lower temperatures. The semi-log plot indicates 
that an exponential form r ~ exp (k^/hsT) 1 *^, with £ = 0.7, describes the data well in the 
entire range. Though such a dependence is expected according to the random first order 
theory (RFOT) [?| , the exponent value we observe cannot be easily rationalized within that 
framework. We comment further on the significance of the exponent value later. Figure 3 
shows relaxation times r(T, N) for different N values scaled to the asymptotic N — ► oo value 
r(T), plotted against values of Xi(T,N) scaled to the asymptotic N — ► oo value x\(T)- If 
the system size dependence of r and x\ are governed by the same length scale, one must 
expect a universal dependence of the scaled relaxation times on the scaled x\ values. From 
the data shown in Figure 3, it is clear that there is no systematic relation between the scaled 
r and x\ that describes their variation both with T and with N. These results indicate that 
the growth of r with decreasing T is not governed solely by the correlation length £ that 
describes the growth of Xi- 
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Motivated by the AG relation 1], r oc exp where A is a constant, we next con- 

sider the dependence of r on the configurational entropy S c whose evaluation is described 
elsewhere 24J. As shown in Figure 4 where log(r) is plotted vs. for all temperatures 
and system sizes studied, we find a remarkable agreement with the AG relation, not only 
vs. T but also for all system sizes. To our knowledge, such a demonstration of the valid- 
ity of the AG relation for finite or confined systems has not been made earlier. Thus, the 
iV-dependence of r, which can not be understood from dynamical finite-size scaling, can 
be explained in terms of the iV-dependence of S c , suggesting that the growth of r with 
decreasing temperature is more intimately related to the change of S c , than to the increase 
of the correlation length £ and susceptibility \ii associated with dynamical heterogeneity. 
As S c at a given temperature varies with system size N, it is tempting to inquire if the 
iV-dependence of S c is associated with a length scale. We extract such a length scale from 
data collapse of S C (T, N), scaled to its value as N — > oo, shown in the top inset of Figure 4. 
We obtain reasonable data collapse, but the extracted length scales turn out to have sub- 
stantially weaker T-dependence compared to £. The comparison of the length scale obtained 
from the configurational entropy and that obtained from dynamical heterogeneity is shown 
in the bottom inset of Figure 4. 

A central role for the configurational entropy, along with an analysis of a length scale 
relevant to structural relaxation, are the content of the random first order theory, developed 
by Wolynes and co-workers Q. According to RFOT, the length scale of dynamical hetero- 
geneity is the "mosaic length" that represents the critical size for entropy driven nucleation 
of a new structure in a liquid. Mean-field arguments based on known properties of infinite- 
range models suggest that the RFOT mechanism is operative for temperatures lower than 
Tmct- In this regime, the dynamics of the system is activated, with the relaxation time 



expected to vary as r = tq exp 



, where AF is the free energy barrier to structural 



rearrangements, and ip is an unknown exponent. The free energy barrier in turn depends on 
the mosaic length as ~ £ 61 , where 9 describes the dependence of the surface energy on 
the size of a region undergoing structural change. The observed validity of the AG relation, 
and the dependence of the relaxation time on the length scale £, r ~ exp (k^/ksT)^, with 
( = 0.7, can be rationalized within RFOT if the exponent 9 is assumed to be close to 2.3, 
and the exponent ip is close to 0.3. However, this interpretation has the drawback that the 
exponent 9 does not satisfy the physical bound, 9 < 2, in three dimensions, and there is no 
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evident explanation for the value of ip. We note that similar conclusions were reached in 



a recent analysis 



2l| of experimental data near the laboratory glass transition, on a large 



class of glass-forming materials. Thus we find puzzling values for the exponents relevant to 



the applicability of RFOT, which are in need of explanation, and data in [21| indicate that 
such a result may apply for a wide range of temperatures, all the way to the experimental 
glass transition. 

RFOT focuses on behavior near the glass transition, and in the limiting case of the spin 
glass models where theoretical perditions are available, configurational entropy plays no 
role in the behavior of the system above the mode coupling temperature. However, there 
have indeed been attempts to extend the RFOT analysis to temperatures above the mode 



coupling temperature 



31 



321 ] and to estimate a mosaic length scale at such temperatures, anc 
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we thus compare our results with predictions arising from these analyses. Stevenson et al 
have considered the change in morphology of rearranging regions above the mode coupling 
temperature, and correspondingly the dependence of relaxation times on configurational 
entropy. The predicted dependence of relaxation times on configurational entropy differs 
from the Adam-Gibbs form, whereas our results strikingly confirm the Adam-Gibbs form. 



Franz and Montanari 



32J have estimated a mosaic length scale in addition to a heterogeneity 



lenght scale, and have discussed the crossover in the dominant lengthscale near the mode 
coupling temperature. The significantly smaller mosaic length scale they estimate appears 
consistent with the smallness of the length scale we obtain from the finite size dependence of 
the configurational entropy. Nevertheless, this analysis does not contain explicit predictions 
regarding the relevance of the configurational entropy at temperatures higher than the mode 
coupling temperature. 

Our observation that the configurational entropy predicts the relaxation times in accor- 
dance with the AG relation for all the temperatures and system sizes we study poses serious 
challenges to current theoretical descriptions based on the analogy with the behaviour of 
mean-field models. Although the relevance of the configurational entropy at high temper- 
atures has been observed in earlier simulation studies and analyses based on the inherent 



structure approach 



24 



33 



34J, we emphasize that a theoretical analysis that satisfactorily 



explains such dependence is not at hand at present, and our results concerning the robust- 
ness of the Adam-Gibbs relation in finite systems highlights further the challenge to existing 
theoretical descriptions. Equally importantly, our results reveal that the length scale asso- 
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ciated with dynamical heterogeneity does not play the central role attributed to it in recent 
analyses, and highlights the necessity to understand the role of other relevant length scales, 



along the lines of the analysis in 32]. 
Methods : 

Simulation details: The system we study is a 80:20 (A:B) binary mixture of particles inter- 
acting via the Lennard- Jones potential: 

V ad r,^ eas [(^) n -^f)"], (1) 

where a,/3e {A, B} and, e A B^AA = 1-5, e BB /e A A = 0.5, oabI^aa = 0.80, cr BB /(T AA = 0.88, 
masses m A = m B . The interaction potential is cutoff at 2.50<t Q/ 3. Length, energy and time are 



reported in units of o~aa, £aa and y o~ 2 AA / 'caa, and other reduced units are derived from these. 
All simulations are done for number density p = 1.20. We have used a cubic simulation box 
with periodic boundary conditions. Simulations are done in the canonical ensemble (NVT) 
using a modified leap-frog integration scheme. We simulate for seven temperatures in the 
range T G {0.450, 1.00}. The mode coupling temperature for this system has been estimated 



23j to be T MC t — 0.435. We equilibrate the system for ~ 10 7 — 10 8 MD steps depending on 
system size and production runs are at least 5 times longer than the equilibration runs. We 
use integration time steps dt from 0.001 to 0.006 for the temperature range 0.800 to 0.450. 
The studied system sizes vary from A^ = 50 to A^ = 1600. 

Dynamics: Dynamics is studied via a two point correlation function, the overlap Q(t), 

N 

Q(t) = I drp(r,t )p(r,t + t )~J2w(\n(t )-ri(to + t)\) (2) 

i=l 

where p(r,t ) etc are space-time dependent particle densities, w(r) — 1, if r < a and zero 
otherwise, and averaging over the initial time to is implied. The use of the window function 
[a = 0.30] treats particle positions separated due to small amplitude vibrational motion as 
the same. The second part of the definition is an approximation which uses only the self- 
term, which we have verified to be reliable (see for details). The structural relaxation 
time r is measured by a stretched exponential fit of the long-time decay of Q(t). 
The fluctuations in Q(t) yields the dynamical susceptibility: 

X4(f) = ^l(Q 2 (t)} - (Q(t)} 2 }. (3) 

Previous work has shown that X^(t) reaches a maximum for times T4 which are propor- 
tional to the structural relaxation time r. We report the values of Xa = Xiit — T i)- 



The Binder cumulant, which we use for finite size scaling, is defined as 

r(nt) < [QfrO- < Q(n) >] 4 > 

B{N ' T) ~ 3 < IQ(t 4 )- < Q(n) >] 2 > 2 " (4) 

B(N,T) = 0, if the distribution P(Q(t4)) is Gaussian, and is a scaling function of £/L 
only (where L is the linear size of the system, and £ is the correlation length), without any 
prefactor. 

Configurational Entropy : S c , the configurational entropy per particle, is calculated as the 
measure of the number of distinct local energy minima, by subtracting from the total entropy 
of the system the 'vibrational' component: 

S c (T) = S total (T)-S mb (T). (5) 



Details of the calculation procedure are as given in 24|. 
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Figure Captions 

Figure 1: System size dependence of dynamic susceptibility x\i an d finite size 
scaling of the Binder cumulant B(N,T). a Inset: x±(f)i shown here for N = 1000 
and selected temperatures, exhibits non-monotonic time dependence, and the time r 4 at 
which it has the maximum value has been observed to be proportional to the structural 
relaxation time r. Main panel: Peak height of the four-point dynamic susceptibility, 
Xi(T,N) = x^it = t a^T,N), has been shown as a function of system size N for different 
temperatures. For each temperature, Xi(T, N) increases with system size, and saturates for 
large system sizes. x\{T, N) also increases as the temperature is lowered, b Inset: The 
distribution P[Q(r 4 )— < Q{t$) >] of Q{j±) — < Q{r±) > is shown for two system sizes for T = 
0.470. While the distribution for the large size is nearly Gaussian, the small system exhibits a 
strongly bimodal distribution. Such bimodality is also observed to emerge as the temperature 
is decreased at fixed system size. Main Panel: Binder cumulant B(N,T) (see Methods ) has 
been plotted as function of iV/£ 3 for different temperatures in the range T e [0.45, 0.80]. The 
correlation length £ is an unknown, temperature dependent, scaling parameter determined by 
requiring data collapse for values at different T. By construction B(N, T) will approach zero 
for large system sizes at high temperatures. It changes to negative values as the temperature 
or the system size is decreased such that -P[<5(r 4 )— < Qij^) >] becomes bimodal. The 
correlation length £ (T) is the only unknown to be determined in order to obtain data collapse 
for B(N,T) and the quality of the data collapse confirms the reliability of this procedure. 

Figure 2: Relaxation times as a function of temperature and system size. 
Relaxation time r(T, N) for the largest system size increases roughly by three decades from 
the highest to the lowest temperature shown. For each temperature, r(T, N) increases as 
./V is decreased for small values of N, displaying a trend that is opposite to that observed 
near second order critical points. For the smallest temperature, r(T, N) increases by about 
a decade from the largest to the smallest system size. 

Figure 3: Relationship between the relaxation time r(T, N), correlation length 
£(T) and the dynamic susceptibility Xa{T, N). Inset: (a) t(T,N — > 00) is shown 
against £(T), in a log-log plot (bottom curve). This plot shows that a power-law dependence 
holds over a temperature range above T = 0.5, but breaks down at lower temperatures, (b) 
t(T,N — > 00) is shown against £(T)//csT, in a semi-log plot (top curve). This plot shows 
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that an exponential dependence r ~ exp ^/c(^//c B T)^j, with ( = 0.7, describes the data well 
in the entire temperature range. However the observed exponent value ( = 0.7 is difficult 
to explain with existing theories. Main Panel: Relaxation times t(T,N)/t(T,N — > oo) 
shown against ^(T , N) / ^(T , N — > oo) in a semi-log plot. While at fixed N both r and 
X4 increase upon decreasing T, at fixed T, they show opposite trends, with r increasing for 
decreasing N and \\ increasing for increasing N. If r and x*l are determined by the same 
length scale £ and further, if their finite size behavior is governed by iV/£ 3 , the plotted data 
are expected to lie on a universal curve, which is seen not to be the case. 

Figure 4: The dependence of relaxation times r on the configurational entropy 
Sc. Main Panel: The relaxation times are shown against configurational entropy in an 
"Adam-Gibbs" plot [log(r) vs 1/(TS C )}, for all the temperatures and system sizes studied. 
The impressive data collapse of all the data onto a master curve indicates that the configu- 
rational entropy is crucial for determining the relaxation times. The overall behavior is well 
described by the AG relation, which requires log(r) oc 1/(TS C ). Although other powers of 
1/(TS C ) may also describe the data well, the improvement in fit quality is marginal, and 
hence we treat the data presented as validating the AG relation. Top Inset: The configu- 
rational entropy S C (T, N) scaled to its N — > oo value, has been plotted as function of N/^ 
for different temperatures in the range T e [0.45, 0.80], to extract a temperature dependent 
length scale £ s that leads to data collapse. Bottom Inset: The length scale obtained from 
the data collapse of the configurational entropy (green diamonds) compared with the length 
scale obtained from finite size scaling of the Binder cumulant (red triangles). It is apparent 
that the length scale from configurational entropy shows very weak temperature dependence, 
in contrast with the dynamical heterogeneity length scale. 
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